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^ ■ Synopsis 

^ ! We consider baryon vertices within the gauge/gravity correspondence for a class of curved 

p 

Q ■ allows classical solutions representing bound states of /c-quarks with k less than or equal 



OO 
O 



X 



to A^. We construct the corresponding classical configurations and perform a stability 
analysis. We present the details for the theory at the conformal point and at finite 
temperature and show that there is a critical value of k, below which there is instability. 
This may also arise when the baryon reaches a critical size. We also extend our treatment 
to magnetically charged baryon vertices. 



Contents 

1 Introduction 1 

2 The classical solution 3 

3 Stability analysis 7 

3.1 Small fluctuations 7 

3.2 Boundary and matching conditions 9 

3.3 Zero modes 10 

4 Examples 11 

4.1 The conformal case 11 

4.1.1 Classical Solution 11 

4.1.2 Stability analysis 12 

4.2 Non-extremal D3-branes 13 

4.2.1 Classical solution 13 

4.2.2 Stability analysis 14 

5 Magnetic flux and stability of the baryon vertex 15 

5.1 The conformal case 16 

5.2 Non-extremal DS-hianes 17 

6 Summary and Discussion 18 

1 Introduction 

The basic construction in use to understand baryons in the context of the gauge/gravity 
correspondence [1] involves A^ heavy external quarks on the boundary of AdS^ at the end 
points of A^ fundamental strings terminating at a vertex on the interior of AdS^. Each 
string contributes a unit of charge, so that for the total charge to be conserved at the 
vertex, it was proposed in [2] that a D5-brane wrapped on the internal S'^-sphere should 
be placed there. Due to its presence, a coupling via a Wess-Zumino term to the Dirac- 
Born-lnfeld (DBl) electric field strength ensures that the total charge indeed is preserved. 



This idea was actually applied in calculating the classical solution corresponding to such a 
configuration in [3, 4] and involved generalizations of techniques developed for the heavy 
quark- ant iquark system at the conformal point of A/" = 4 SYM for SU{N) in [5]. 

In recent works [6, 7, 8] we have extensively studied within the gauge/gravity correspon- 
dence, the stability of interaction potentials between a heavy quark and anti-quark in 
mesons, as well as between a quark and monopole, for a general class of backgrounds. A 
general outcome of these works is that a stability analysis is very important in establish- 
ing the physical parameter space of the theory which is typically reduced in comparison 
with that allowed from simply the requirement that the classical solutions exist. A sec- 
ond lesson is that in cases of multibranch solutions the point where instability begins 
does not always coincide with the point where a solution bifurcates. It seems that this 
feature requires the coexistence in the classical configuration of fundamental as well as of 
solitonic objects. In particular, this was the case of the potential for the quark-monopole 
pair in which a string junction was used that requires precisely a fundamental-, a D- and 
a dyon-string. 

The case of baryons we focus in this paper presents a new challenge and requires special 
care. In particular, at a technical level it requires a multistring configuration with A^ 
strings. In addition, the vertex where these meet in the bulk of AdS^ is not a mathemat- 
ical point, as in the case of a string junction, but it is a place where a solitonic object, 
the D5-brane wrapped on the internal S^, is placed. A further challenge, that initially 
prompted our analysis, is that there have been configurations with 5A^/8 < k < N quarks 
[3, 4] that are perfect as classical solutions to the equations of motion, besides the ex- 
pected ones with k = N. This fact is against physical intuition and experimental results 
that, bound states of quarks and antiquarks are singlets of the gauge group. Of course in 
the case at hand we deal with a maximally supersymmetric theory, but one would like to 
show that a better bound for k/N can be found, if not k = N exactly. It is the purpose 
of the present paper to address these and related issues. 

The organization of this paper is as follows: In section 2 we formulate baryons for a 
general class of backgrounds. In particular, we derive general formulas for the baryon 
binding energy for the case of quarks uniformly distributed on a spherical shell when 
forming the baryon. In section 3 we perform a stability analysis under small fiuctuations 
and derive their equations of motion as well as the boundary and matching conditions 
at the boundary of AdS^ and at the baryon vertex. In section 4 we present examples of 



baryon configurations within the gauge/gravity correspondence, using the conformal and 
finite temperature backgrounds. We will show that stability restricts the value of k to be 
larger than a critical value even in the conformal case (a higher than 5A^/8 value). Also 
even a part of the energetically favorable branch of multi-branched potentials should be 
disregarded as unstable. In section 5 we extend our treatment to the case of magnetically 
charged baryon vertices and show that stability requires un upper bound for the instanton 
number, as compared to the 't Hooft coupling, associated with a self-dual gauge field. 
Finally, in section 6, we summarize our results and discuss some future directions. 

2 The classical solution 

In this section we develop and present the general setup of the gauge/gravity calculation 
of the potential energy of a baryon consisting of heavy static quarks. Since the underlying 
dual theory is A/" = 4 SYM for SU{N), these computations involve A^ strings as well a 
D5-brane wrapped on the internal S'^-sphere. Certain elements of the computation will 
be based on the well known techniques used for the calculation of the potential of a heavy 
quark- ant iquark pair [5] as generalized beyond the conformal point in [9], as well as the 
concept of the baryon vertex [2]. Our classical construction will essentially generalize 
that for the conformal case in [3, 4] to a more general class of backgrounds. 

We consider diagonal metrics of Lorentzian signature of the form 

rfs2 ^ Q^^^^2 ^ G^,^{dx^ + dy"^ + dz^) + Guudu^ + R^d^l , (2.1) 

where x, y and z denote cyclic coordinates and u denotes the radial direction playing the 
role of an energy scale in the dual gauge theory. It extends from the UV at -u ^ oo down 
to the IR at some minimum value Mmin determined by the geometry. 

It is convenient to introduce the functions 

j\u) = —GttGxx I dW = ~GttGuu ) h[u) = GyyGuu ■ (2-2) 

In the conformal limit we can approximate the metric by that for AdS^ x S^ with radii 
{a' = 1) R = (AjcgsNy^*, with gg being the string coupling. In this limit the leading 
order expressions are 

f{u) ~ u^ , g{u) ~ 1 , h{u) ~ 1 . (2.3) 
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Figure 1: A baryon configuration with k-exteinal quarks placed on a spherical shell of 
radius L at the boundary of AdS space, each connected to the D5-brane which is wrapped 
on an S^ located at u = Uq and N — k straight strings ending at Mmin- 

A baryon is by definition a bound state of A^-quarks which form the completely antisym- 
metric representation of SU{N). However, within the gauge/gravity correspondence, we 
can in general construct a bound state of fc-quarks with k < N. This bound state consists 
of a D5-brane located in the bulk and fc-external quarks at the boundary of AdS^. Each 
quark is connected with a string to the DS-brane. There are also N — k straight strings 
that go from the D5-brane straight up at u^m- This description is depicted in Fig. 1. 
We will take N,k ^ 1, such that the ratio k/N is held finite. 

At this point we have to choose a profile for the quark distribution inside the baryon. 
We expect a radially symmetric distribution and for simplicity we take it to be a uniform 
distribution on a spherical shell of radius L. We will briefly comment on other choices 
and on how the preferable distribution could be determined in the last section. Also 
note the entire description is valid only if the endpoints of the A^-strings are uniformly 
distributed on the internal S^-sphere so that the latter is not deformed and the probe 
approximation still holds (for a discussion on this point see [4]). Such a choice completely 
breaks supersymmetry [10, 11]. 

Within the gauge/gravity correspondence, the binding potential energy of the baryon is 
given by 



g-iET _ gi5., 



(2.4) 



where Sd is the classical action of the above configuration. This action consists of three 
terms, a Nambu-Goto action for each of the A^-strings, a Dirac-Born-Infeld action for 
the D5-brane and a Wess-Zumino term between the A^-strings and the DS-brane 



Sng[C] = tt drda^-deigc^p , 

^T^ Jc 



^D5 = 7T-^r- / d^V-h^J - 27rF,, , (2.5) 

Swz = , !, I d^x C(^) A F , M = RxS^ , 
{2'nygs Jm 

where gafs-, hij, Fij are the induced metrics for the string, D5-brane and Born-Infeld field 

strength, respectively. We will consider the simplest case where the DBI-action is just 

given by its induced metric and the presence of the Wess-Zumino term acts as a charge 

conservation term [2]. However, there are cases of non-trivial Fij [12] to which we will 

return later in section 5. We have not also included the contribution of the antisymmetric 

NS-NS two-form due to the fact that it is not present in our examples. 

We first fix reparametrization invariance for each string by choosing 

t = T , u = a . (2.6) 

For static solutions we consider the embedding of the S'^-sphere on the D3-brane in 
spherical coordinates (r, 9, 0) 

r = r(u) , (j), 9 = const. , S'^-angles = const. , (2.7) 

supplemented by the boundary condition 

u{L) = oo. (2.8) 

Then, the Nambu-Goto action reads 

S = -^jdu^Jg{u) + f{u)r'\ (2.9) 

where T denotes time and the prime denotes a derivative with respect to u. From the 
Euler-Lagrange equations of motion we obtain 

^'^^ -fl'' =^ r'd = ^, (2.10) 



where Ui is the value of u at the turning point of each string, /i = f(ui), /o = /(mq) and 

F = Y^^- (2.11) 
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The N — k straight strings which extend from the baryon vertex to -u = -Umin are straight, 
since r = const, is always a solution of the equations of motion (with /i = 0). In the 
examples to follow Umin ^ Ui ^ uq and f{u) is an increasing function. Integrating (2.10), 
we can express the radius of the spherical shell as 



7i / du^ . (2.12) 



MO 



Next we fix the reparametrization invariance for the wrapped D5-brane by choosing 

t = T, 9a = a^, a = 1,2,..., 5. (2.13) 



For the ansatz given above, the DBI action reads 

TNR 



Sd5 — —^ — y —Gtt 

OTT u=uo 

where we have also used that the volume of the unit-S*^ is tt^. 



(2.14) 



27r 1 


/ duVF- 


/ du^/g + 

k 


1 


-a r ^ ^ , R^-Gu 
a X.'"^+ 4a 

< a ^ 1 . 



Finally, inserting the solution for r^j into (2.9), subtracting the divergent self-energy 
contribution of disconnected worldsheets and adding the energy of the D5-brane (2.14), 
we write the binding energy of the baryon as 

k { f^ f^ 1 — — 

E= — {l duVf- I du^+^^ I du^+ ^'^ ,y ), (2.15) 

where 

(2.16) 

The expressions for the length and the energy, (2.12) and (2.15) depend on the arbitrary 
parameter ui which should be expressed in terms of the baryon vertex position uq. The 
most convenient way to find this is to impose that the net force at the baryon vertex 
is zero. Following an analogous procedure to that for the string junction in [8], we 
temporarily change gauge to r = a and from the boundary terms arising in deriving the 
classical equations of motion we find that for a uniform distribution of the quarks on a 
spherical shell, we get the relation (reducing to that in [3] for the conformal case) 

1 -a R 



COS0 = h - — -du\/-Gtt 

a ^o-^g u=uo 

cose = v/l-/i//o, (2.17) 

where O is the angle between each of the /c-strings and the -u-axis at the baryon vertex 
and which determines Mi in terms of uq. The above relation can also be interpreted as 
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the no-force condition for its -u-component [4]. The term on the left hand side of (2.17) 
is due to the strings attached to the boundary at the UV. The force they exert on the 
baryon vertex is balanced by the forces exerted on that by the straight strings and the 
tendency of the Z^S-brane to move towards u^[^, represented by the first term and the 
second term on the right hand side of (2.17) (the second term is indeed positive in our 
examples). Hence, if a becomes small enough we expect that this condition might not be 
possible to satisfy. As we will see in the examples to follow, (2.17) has a solution for a 
parametric region of (a,-Uo)- However, in order to isolate parametric regions of physical 
interest a stability analysis of the classical solution should be performed, which, as we 
will see, further restricts the allowed region. 

3 Stability analysis 

We now turn to the stability analysis of the above classical configuration, aiming at 
isolating the physical parametric regions. This involves string fluctuations which fall into 
the general results analyzed in [6, 7] for a single string. However, the quantitative results 
are quiet different, due to the relative orientation of the strings, being distributed on a 
sphere, and the different boundary conditions that will be imposed here. 

3.1 Small fluctuations 

We parametrize the string fluctuations about the classical solution, by perturbing the 
embedding according to 

r = rci{u) +6r{t,u) , 9 = Oq + 69{t,u) , (p = (po + S(f){t,u) . (3.1) 

We have kept the gauge choice (2.6) unperturbed by using worldsheet reparametrization 
invariance. We then calculate the Nambu-Goto action for this ansatz and we expand it 
in powers of the fluctuations. The resulting expansion for the quadratic fluctuations for 
each one of the /c-strings can be written as 



5*2 = / dtdu 

271 



nf h i'r"^ hr"^ F^/"^ 



2F3/2 2F1/2 2F1/2 2g 

frl sinHo ^,,2 hrlsin^oF^/^ ■ 



2F1/2 ^ 2g 

Using the Euler-Lagrange equations and the ansatz 



(3.2) 



(5x^(t, u) = 6x^'{u)e~'''' , x^ = r, ^, , (3.3) 



we find for eacli of the /c-strings tlie linearized equations for tlie longitudinal and trans- 
verse fluctuations 
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du 


( 9f d\ 

[ry^duj 
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' d 
du 


( frl d \ 
\F^l^duJ 
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' d 
du 


{ hi d \ 
\F^l^duJ 
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(3.4) 



For the straight strings the expressions for the action and the equations of motion were 
given in [8] and will not be repeated here. As discussed extensively in [6, 8] instabilities 
are indicated by the presence of a normalizable zero-mode which can exist for speciflc 
values of uq, denoted by uqc- Such values are in general different than possible values 
of Mo at which the length reaches an extremum, denoted by -uom and found by solving 
L'{uo) = 0. Our aim is to flnd for the zero mode the critical curve in the parametric 
space of a and Uq which separates the stable from the unstable regions. 

Finally, we turn to the fluctuations of the D5-brane and perturb the embedding according 
to 



x'^ = 5x^{t, 9a) , U = Uq , X^ = X, y, z 



(3.5) 



leaving the position of the D5-brane at u = Uq intact due to the gauge choice u = a for 
the strings. The second order, in the fluctuations, term in the expansion of the D5-brane 
action is 



S 



NR 



57r" 



G,, 



G. 



dtdVt^^^-Gtt{ 1 + -^^I'^^dJx^'dpSx'' + -^Sx^'dx" 

2 2Gii 



(3.6) 



where 7q,/3 is the metric of S^ and the action is calculated at -u = uq. The subscripts 
Q;,/i refer to the angles of S^ and to the x,y,z coordinates, respectively. We expand 
the fluctuations in terms of the spherical harmonics of the S'^-sphere "$£, satisfying the 



eigenvalue eq. V^\E'^ 



4)^,, £ = 0,1,2,..., as 
6x^'{t,9a) =6x^'{t)^e{9a) 



Then, from the Euler-Lagrange equations for the action we flnd that 

+ fi^2(5x^ = , nj = -GttiuoW + 4) > . 



dt^ 



(3.7) 



(3.^ 



Thus the classical solution is stable under the 6x'^ fluctuations. Note that, there are 
no boundary conditions for these fluctuations, the reason being that the M x S*^ has no 
boundary. 



3.2 Boundary and matching conditions 

To fully specify our eigenvalue problem, we must impose boundary conditions on the 
string fluctuations in the UV limit, -u ^ oo and in the IR limit, u = u„^i^, as well as 
matching conditions at the baryon vertex, i.e. u = Uq. The boundary condition at the 
UV are chosen so that the quarks are kept flxed, that is 

$(«) = , as M ^ oo , (3.9) 

where $ is any of the fluctuations of any of the /c-strings. At -u = -Umin we demand 
finiteness of the solution and its -u-derivative for the N — k straight strings that extend in 
there, so that the perturbative approach is valid. It turns out that for the zero-mode they 
completely decouple from the analysis, as in the case of the quark-monopole potential in 
[8]. Since we consider a uniform distribution of quarks on a spherical shell, the D5-brane 
will be located by symmetry aX x = y = z = 0. For each one of the /c-strings, 6r is its 
longitudinal fluctuation and 66, 6(f), are its transverse ones. 

For the longitudinal fluctuations 6r, it is more convenient to work out the necessary 
boundary condition in the gauge r = a, instead of -u = o" and then translate the result in 
the latter gauge where the equations of motion are simpler. In that spirit, following also 
an analogous treatment in [6, 8], the gauge change is effectively done by the coordinate 
transformation 

u = u + 6u{t,u) , 6u{t,u) = — !1-1^ . (3.10) 

Hence, the boundary condition at the baryon vertex can be found, if we temporarily 
change gauge to r = o" and apply small fluctuations for the w-coordinate. Demanding a 
well behaved boundary problem for (5-u-fluctuations we obtain 

d6u 



dr 



= 0. (3.11) 

From (2.10), (3.10) and (3.11) we flnd that the (5r-fluctuations, for each of the fc-strings, 
satisfy the following boundary condition 

u = Uo: 2(/ - h)6r' + 6r Uf - jrfi - ^(/ - fi)) = . (3.12) 

This is an equation for Uq in terms of a and the parameters of the speciflc problem. 
We will denote its solution, whenever it exists, by -uqc- We will not assume continuity 
relations between the fluctuations of the D5-brane and the A^-strings. 
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In general, the value of mqc does not coincide with that corresponding to a possible 
extremum of the radius L{uo), denoted by -Uom- For the quark-antiquark potential these 
values coincide [6, 7], but they do not do so for the case of the quark-monopole potential 
[8]. In our case using (2.12), (2.15) and (2.17) we find that 

dE _ k^/hdL 
duo 271 duQ 



where from (2.17) we compute that 



^^9^ = ""®^°-2^^°^°^®^H^^^J 



(3.14) 



A necessary condition for -uom to exist is that the right hand side of (3.14) is positive, 
which is indeed satisfied in the non-extremal D3-branes example we will examine in the 
next section. Also the extrema of L{uq) and E{uq) at uq = uom coincide as in the case of 
the quark-antiquark and quark-monopole potential energy interactions shown in [7] and 
[8], respectively. Similarly to [9], we can easily show that 

dE__^nr 
'dL~ 2^^^^ ' 

fE_ ^ J^Ji 1_ dui 

dL^ 47r yfj; L'{uo) duo' 

Since dE/dL > the force is manifestly attractive. On the other hand its magnitude 
may decrease or increase with distance depending on the sign of the right had side of the 
expression for d'^E/dL'^. Note that this is not necessarily in conflict with gauge theory 
expectations since there is no analog for the binding energy of baryons of the concavity 
condition of the heavy quark-antiquark pair found in [18]. In addition, having assumed 
that the quarks, the baryon consists of, are located on a spherical shell has an inherit 
degree of arbitrariness that is reflected in the details of the expressions we have derived. 

3.3 Zero modes 

We will now study the zero mode problem for a spherical shell distribution and flnd the 
critical curve in the parametric region of (a, -Uq)- 

The longitudinal and transverse fluctuations of each string, satisfying already the condi- 
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tions in the UV, are 

5r = A J{u) , 59 = B K{u) , 




where A,B and C are integration constants. Using (2.3) and (2.10) we see that K{u) 
diverges, thus we should choose B = C = 0. Therefore, the transverse fluctuations are 
stable as in the cases of the quark- ant iquark and quark-monopole potentials [6, 8]. 

For the longitudinal fluctuations we should substitute the above expression for 6r into 
(3.12) and obtain a transcendental equation that defines the critical curve separating the 
stable from the unstable regions. 

4 Examples 

In this section we will first review the behavior of the baryon potential emerging in 
the conformal case [3, 4] and examine its stability behavior. Then we will present the 
analogous analysis for the finite temperature case. 

4.1 The conformal case 

4.1.1 Classical Solution 

At first we will consider the conformal case, in which the metric of AdS^ x S^ has the 
form of (2.1), with 

—Gtt = Gxx = Guu = "^2 ■ ^^-^^ 

The radius and the energy corresponding to a uniform distribution of quarks on a spher- 
ical shell are given in terms of the position of the D5-brane uq and the turning point Ui 
of each string 

R'^ul ^ ^ _ kuo ( ^ 5 - 4a' 



^ = 1^^- ^ = i7-^+^;ri' <^-^' 



with 



^-^iM-H)- ^-^(-H'fii ("' 
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and where 2-^1(0, b, c; z) is the hypergeometric function and uq^ Ui^ 0. From (2.17) we 
find that in this case the no-force condition on the w-axis yields 

u, = Uo{l-\'f'\ A = ^. (4.4) 

4a 

We note that (4.4) does not restrict the domain of -uq? but this will not be the case 
when temperature is turned on. Since A < 1, a baryon configuration exists if a > a<, 
a< = 5/8 = 0.625. This value was found in [3, 4]. Using (4.2) and (4.4) we can find the 
binding energy in terms of the physical length 



^_^<.Vr^/ 5-4aA 

27rL 3 V 4a y ^ ^ 

From this expression it is apparent why the possibility a = 5/8 (A = 1) is discarded. The 
energy of the baryon goes as 1/L as it was expected by conformal invariance and has the 
familiar for these type of computations dependence on the 't Hooft coupling. 

4.1.2 Stability analysis 

Next we consider the stability behavior of the conformal case. For the longitudinal 
fiuctuations we have to solve (3.12), where 

Substituting (4.4) and (4.6) in (3.12) we find the following transcendental equation for 
the zero mode to exist 

Using (4.4) and (4.7) we find numerically a critical value for a, i.e. Qc — 0.813 . 

Thus, the stability analysis sets a low bound for a which is still less than unity. However, 
given that this corresponds to the conformal point of a maximally A/" = 4 supersymmetric 
theory we find it remarkable. The only other case of perturbative instability in the 
conformal point is for the quark- ant iquark potential in the /9-deformed A/" = 1 theory 
[16], whose supergravity dual was constructed in [17]. In this case the instability occurs 
for values of the cr- deformation parameter larger than a critical value [7] . 
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4.2 Non— extremal D3— branes 

4.2.1 Classical solution 

In this section we consider a stack of A^ non-extremal D3-branes [13]. In the field theory 
limit the metric has the form of (2.1) where 

Zi ( u \ u 

<^« = -G'«« = -;^ ( 1 - :^ ) 1 Gxx = ^- (4.8) 

The corresponding Hawking temperature is Th = ^/{tiR^). In what follows, we will 
compute in units of /i, so that we set /i = 1. The radius and the energy read 






1 1 1 3_ 1 wf 

~4' 4 

where -Fi(a, 6i, 62, c; Zi, Z'l) is the Appell hypergeometric function and -Uq ^ "^i ^ 1- From 
(2.17) we find that in this case the no-force condition on the -u-axis yields 



^ = -"°^^'-4'-2'2'4'^'^'+' 



^ 1-a 1 + I/mq ^ \uq-u\ ,,^r.\ 

cos0 = \ , ' " cose=W-5i^ i. (4.10) 

Employing the fact that cos© ^ 1 in (4.10) we end up with a second order polynomial 
inequality which has a solution for a ^ 5/8, namely 



Wo ^ 'U<(a) , ■u<(a) = 1/A 



1/4 



A = -9 - 32a(a - 1) + 4V2(2a - 1)^3 - 8a(l - a) , (4.11) 

a G [5/8,1], A G [0, -9 + 4^6] , mq e [00, «<(!)] , m<(1) ~ 1.058 . 

However, if a = 5/8 there is no classical solution since -uq -^ 00. For a > 5/8 the solution 
of (4.10) is 

u, = uo{l - A\uo)Y/' , A(..o) = ^^a/T^Tk + ^^ . (4.12) 

In this case we can not express E in terms of L. In fact i5 is a doubled-valued function 
of L with energetically favorable and unfavorable branches, as depicted in Fig. 2a. It is 
similar to the quark- ant iquark potential for the finite temperature case [14]. 
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Figure 2: (a) Plot of E{L) for non extremal Z^S-branes. The various types of line cor- 
respond to stable (solid grey) and unstable (dashed grey) configurations. As a decreases 
unstable part extends to the lower branch as well. This behavior is described in section 
(4.2.2). (b) From the left to the right: i) Plot of m< versus a, ii) Plot a versus Wom, hi) 
Plot of Uqc versus a. The regions on the right and left hand side of the stability line are 
stable and unstable, respectively. 

4.2.2 Stability analysis 



Next we will consider the stability behavior of the non-extremal D3-branes. For the 
longitudinal fluctuations we have to solve (3.12) where 

[A 



5r{u) = A 



du- 



V' 



u^ 



:Fi 



3 1 3 7 1 ut 



(4.13) 



'{u^-uff/^ 3m3 '\4' 2'2'4'm4'm4^ 

From (3.12) and (4.13) we find that the zero mode condition is given by the transcendental 
equation 

3 OoA _ „,4 



2a/mo - uj 



2u^-ut-l ^ f3 1 3 7 1 wf 



Mn 



1)3/2 



4' 2'2'4'«^'«4 



(4.14) 



From numerical considerations we find that for every value of a larger than the critical 
value found in the conformal case a^ — 0.813 there is a zero mode on the energetically 
favored branch of E{L) i.e. uqc > Wom- Its solution has the expansion 

0.679 0.778 ,^ / 1 \ , , 

a,^0.813 + -^ + -^ + O -^ . (4.15) 

As Oc decreases mqc increases and for Oc ~ 0.813 even the conformal part becomes unsta- 
ble. The results of our stability analysis are summarized in Fig. 2b. 
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Also, for the particular case of a = 1 we present the table: 



Values for a = 1 


Minimum 


Critical 


Maximum 


Uq 


1.0581 


1.4292 


1.3437 


L/R^ 





0.2920 


0.2951 


E/N 


0.0097 


0.0266 


0.0273 



5 Magnetic flux and stability of the baryon vertex 

In this section we will generalize our discussion to the case where we turn on, besides 
the electric ones, non-trivial magnetic flux components in Fij, as in [12]. The fact that 
this is possible could be seen, from the fact that the S'^-sphere can be regarded as a U{1) 
flbre over CP^ with a non-trivial flbre connection [19]. This f/(l) connection, introduces 
a generalization of the baryon vertex in which non-zero magnetic components of the 
Born-Infeld field F^j appear in the CP^ directions in such a way that it is self-dual with 
respect to that metric. For topological reasons, its flux over non-trivial two-cycles of CP^ 
is quantized and n^ represents the instanton number. Working along similar lines as in 
section (2.2) of [12], we flnd that the action of the D5-brane (2.5) in our conventions gets 
modifled to 

TNR^f^GTt 



S 



D5 



Stt 



where 



u = l + 3271% , b = 



n 
R^ 



(5.1) 



(5.2) 



Note that the relevant parameter is not the instanton number n^ but its ratio with the 
't Hooft coupling. The analogous expression to (2.15) for the binding energy of the 
baryon is 

A; I /•"- , ^ /•"",_ 1 

2^ 



E 



duVF - 



duy/g + 



duy/g + V 



4a 



(5.3) 



U=Uq 



while the expression for the radius of the baryon is given by the same equation as in (2.12). 
The analogous expression to (2.17) for the no-force condition on the w-axis yields^ 

1 



COS0 
COS© 



a , j-,du\/—Gtt 
V u R 



Aa^ 



a/1 - /i//o . 



(5.4) 



^The conditions (3.13) and (3.15) are also valid in this case as well. In using (3.14) one should include 
an overall factor of v to the second term in the right hand side. 
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The stability analysis of this classical solution is easy to discuss. Since the Born-Infeld 
field is self-dual, the determinant under the square root in the action of the D5-brane 
(2.5) is a perfect square, and there no mixing terms of the induced metric of the Db- 
brane and the Born-Infeld field. Furthermore, the Born-Infeld field is a two-form on 
CP that remains unchanged since perturbing the position of the D5-brane we keep the 
gauge choice (2.13). Hence, the stability analysis described in section 3, is applicable for 
non-zero magnetic fiux as well. 

We will next specialize to the conformal and non-extremal D3-branes cases. 

5.1 The conformal case 

We will first consider the classical solution for the conformal case. The whole treatment 
can be done in terms of an effective parameter a^s (defined below) and using our previous 
results in the absence of magnetic charge. The analogous expression to (4.4) reads 

Ml = Mo(l - Aeff)^^^ , Acff = — — , aeflf = 2^ . (5.5) 

4acff 1 + ^i^ 

Since Acg < 1 a baryon configuration exists for 

5 
fleff > 5/8 =^ a > a^ , a< = — h An'^b . (5.6) 

8 

However, since a ^ 1 we get a bound for the instanton number, namely that 

a< < 1 =^ — < ~ 0.0095 , (5.7) 

which is the bound found in [12] for a classical configuration to exist. The binding energy 
in terms of the radius reads 



^^jl/cyr^A^/ _5-4^ (5i 



2txL 3 \ Attcs 

where X and J were given in (4.3). We have shown in section 4 that there is a zero mode 
for acff ^ 0.813, for which the system becomes unstable. This improves the above bound 
for the instanton number, in comparison to the 't Hooft coupling, to 

n 

— < 0.0036 , (5.9) 

that should be respected for the classical configuration not only to exist, but also to be 
perturbatively stable. Stability analysis sets lower and higher bounds for the instanton 
number and a, 0.0036 and 0.813 respectively. This behavior is depicted in Fig. 3. 
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Critical 




Classical 



Figure 3: Classical and critical lines of a versus b, i.e. Qes = 5/8, 0.813, respectively. 
There is a classical solution for values of {a,b), on the right hand side of the classical 
line. As for stability, the regions on the left and the right hand sides of the critical line 
are stable and unstable, respectively. 

5.2 Non-extremal D3-branes 

Finally, we will consider the non-extremal D3-branes. The analogous expression to (4.9) 
for the energy read 



^ k \ ^ I — a , . Un I 1 

E= — {8 + K - 1 + ^T^Wl - - 

ivr 1 a 4a y Mq 

/ 1 1 1 3 1 Mf\ , , 

A + u 

The analogous expression to (4.10) no- force condition, has a solution for a > , 

8 



namely that 



I — a 



1 - 1/u^o + ^ 



1 + iK 



4a 



(5.11) 






- (-8 - z/^ + 32(-a + l)a + 4y2(2a - \)^1-%a{\-a) + z/2 \ . 



The energy versus radius is a doubled valued function with an energetically favorable 
and an energetically unfavorable branch, as depicted in Fig. 2a. 
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Turning now to the stability behavior of the classical solution, we found that for every 
value of a and instanton number above the critical line of Fig.3, there is always a zero- 
mode on the energetically favorable branch of E{L), giving rise to a critical curve bounded 
from below by the conformal critical line. As a and/or the instanton number reach the 
critical line, even the conformal part of the diagram becomes unstable. 

6 Summary and Discussion 

In this paper we first constructed classical configurations with strings and a D5-brane 
that represent bound, baryon-like states, within the gauge/gravity correspondence and 
for a general class of backgrounds. We worked out explicit examples in the A/" = 4 
SYM for SU{N) at the conformal point and at finite temperature. Stability restricts 
the value of k to be larger than a critical value even in the conformal case (0.813A^, 
higher than the classical lower bound 0.625A^). It addition, a part of the energetically 
favorable branch of the multi-branch potential arising in the finite temperature case 
should be disregarded as unstable, even if A; = A^. If the baryon vertex is supplemented 
by magnetic charge, we showed that stability requires un upper bound for the instanton 
number, as compared to the 't Hooft coupling, associated with a self-dual gauge field. 
This is unlike the findings for quark-antiquark potentials in [6, 7] where stability discards 
energetically unfavorable branches completely and for the dyon-dyon potentials in [8], 
where parts of the energetically unfavorable branch become perturbatively stable. 

One may well wonder how generic is our discussion in the finite temperature case, hav- 
ing used a specific black D3-brane background. We have checked that for the Rindler 
space, that represents the horizon of any black hole-type solution similar phenomena 
occur. There is a critical value of k below which there is instability and also part of the 
energetically favorable branch is unstable. Hence, our results are quite generic for a large 
class of backgrounds with black hole-type behavior. Moreover, our techniques could be 
extendable to theories with dynamical flavor, such as the one in [20]. 

Our work is also straightforwardly extendable to supergravity backgrounds in which off- 
diagonal metric elements appear. This cover cases with non-vanishing Wess-Zumino 
term and cases of moving baryons in hot quark-gluon plasmas. The classical solution of 
a moving baryon in hot plasma was studied in [15]. The end result is that for a uniform 
circle distribution at the boundary of AdS^, the binding energy versus length has two 
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branches and for high enough velocities Lmax ~ {'^ — v'^Y^^/Th, where the proportionahty 
coefficient depends on the direction that the baryon moves, with respect to the axis of the 
circle. We have checked the stability behavior and found that for perpendicular to the 
circle distribution motion, there is an instability on the energetically favored branch and 
as speed increases and/or a decreases, even the conformal part of the diagram becomes 
unstable. We expect this behavior to be valid for arbitrary quark distributions in baryons 
and moving directions. It is quite interesting if this is further pursued. 

One could wonder if there are cases when any number of quarks below k = N gives 
rise to unstable configurations, thus realizing our original motivation that prompted this 
work, i.e. only colorless states can exist. It is not a surprise that this does not seem to 
happen for the A/" = 4 SYM theory, since after all this is a maximally supersymmetric 
finite theory. It will be interesting to explore this issue in other theories with reduced 
supersymmetry. 

In our investigation we have chosen to work for simplicity with a quark distribution in 
the baryon that is homogeneous and confined on a spherical shell. This choice is quite 
arbitrary and although it preserves radial symmetry it is not a priory the one that the 
quarks would prefer if, for instance, we apply a variational principle that minimizes the 
total energy of the configuration. It would be very interesting to find the mechanism that 
determines from first principles the appropriate quark distribution in a baryon within this 
framework. As a last comment on that, we have checked that a randomly chosen radially 
symmetric distribution of quarks in the baryon is inappropriate since in that case, not 
even the classical configuration exists, i.e. the no-force condition at the baryon vertex is 
not satisfied. 
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